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ABSTRACT 

We analytically work out the long-term, i.e. averaged over one orbital revolution, time 

variations ( Y \ of some direct observable quantities Y induced by classical and general 

relativistic dynamical perturbations of the two-body pointlike Newtonian acceleration 
in the case of transiting exoplanets moving along elliptic orbits. More specifically, the 
observables Y with which we deal are the transit duration Atd, the radial velocity V p , 
the time interval At cc \ between primary and secondary eclipses, and the time interval 
Ptr between successive primary transits. The dynamical effects considered are the cen- 
trifugal oblateness of both the star and the planet, their tidal bulges mutually raised 
on each other, a distant third body X, and general relativity (both Schwarzschild 
and Lense-Thirring). We take into account the effects due to the perturbations of all 
the Keplerian orbital elements involved in a consistent and uniform way. First, we 
explicitly compute their instantaneous time variations due to the dynamical effects 
considered and substitute them in the general expression for the instantaneous change 
of Y; then, we take the overall average over one orbital revolution of the so-obtained 
instantaneous rate Y(t) specialized to the perturbations considered. In contrast, previ- 
ous published works have often employed somewhat "hybrid" expressions, in which the 
secular precession of, typically, the periastron only is straightforwardly inserted into 
instantaneous formulas. The transit duration is affected neither by the general rela- 
tivistic Schwarzschild-type nor by the classical tidal effects, while the bodies'centrifugal 
oblatenesses, a distant third body X and the general relativistic Lense-Thirring-type 
perturbations induce non- vanishing long-term, harmonic effects on Atd also for circu- 
lar orbits. For exact edge-on configurations they vanish. Both V p and At cc \ experience 
non-vanishing long-term, harmonic variations, caused by all the perturbations consid- 
ered, only for non-circular orbits. Also Ptr is affected by all of them with long-term 
signatures, but they do not vanish for circular orbits. Numerical evaluations of the 
obtained results are given for a typical star-planet scenario and compared with the 
expected observational accuracies over a time span r = 10 yr long. Also graphical 
investigations of the dependence of the effects considered on the semi-major axis and 
the eccentricity are presented. Our results are, in principle, valid also for other astro- 
nomical scenarios. They may allow, e.g., for designing various tests of gravitational 
theories with natural and artificial bodies in our solar system. 

Key words: gravitation-planetary systems-stars: rotation-stars: binaries: eclipsing 



1 INTRODUCTION orbital revolution, precession of the periastron uj due to some 

dynamical perturbations are straightforwardly substituted 
In several papers somewhat hybrid expressions lor certain . . , • , . , , • • , 

in expressions which, instead, contain instantaneous values 
observable quantities pertaining to exoplanets arc presented. r , T , . .,, , . ,„ j. • i i 

° ol the true anomaly f ol the planet. Moreover, a certain lack 

In them, formulas for the long-term, i.e. averaged over one r ... ., . ., ' . , , , , , 

00 ol uniformity in the various treatments tend to make them, 

perhaps, difficult to follow. 

* E-mail: lorenzo.iorio@libcro.it In this P a P er we Propose a more uniform approach 

t Present address: Viale Unita di Italia 68, 70125, Bari (BA), by systematically working out, in an analytic and explicit 

Italy. way whenever possible and/or convenient, the expressions 
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for the long-term variations of various observables induced 
by a number of dynamical perturbations, both of classical 
and general relativistic origin. In particular, we will focus 
on the long-term variations of the transit duration Atd, 
the radial velocity V p , the time interval At ec \ between 
primary and secondary eclipses, and the time Ptr elapsed 
between consecutive primary transits. Concerning gen- 
eral relativity, the issue of the measurability of certain 
eff ects of it in some ext r asolar planets has been tre ated 
in IMiralda-Escudel (I2002T); lAdams fc Laughlinl ll2006allbllch; 




Hevl fc Gladmanl d2007l)~ 



Jordan & Bakos 
Ragozzine fc Wo lf (2009); 



(I2008f): [Pal fc Kocsisl (|20Q, 

iLil (|2010n ; lloriol J2010] ). In particular, the perspectives in 
measuring the corre ction t o the third Kepler law have 



been c onsidered by | lorid (12006 ) and iRagozzine fc Wola 



ll2009h. IMiralda-Escudel (|2002[). faey! fc Gladmanl (|2007h . 
Jordan fc Bakod j20oi) . |Pal fc Kocsisl (|2008h and 
Ragozzine fc Wola (|2009l ) dealt with the possibility of 
detecting the gravitoelectric periastron precessions, while 



ILil (J2010J) looked at the gravitoelectric secular change of 
the mean anomaly A4 conne cted to the v ariation of the 
periastron time transit t p . Adams & Laughlin (2006a be) 
studied the impact of the general relativistic gravitoelectric 
terms on the long-t erm, s e cular interactions among multiple 
planetary systems, lloriol (|20101 ) looked at the relativistic 
effects induced by the rotation of the hosting star. 

In our calculation, we consistently take into account 
the perturbations of all the Keplerian orbital elements in 
the following way. Let us assume that an explicit expression 
is available for a given observable Y in such a way that it is 
function of all or some Keplerian orbital elements, i.e. Y = 
Y(f, {k}), where k denotes the ensemble of the Keplerian 
orbital elements explicitly entering Y apart from the mean 
anomaly A4. Then, we straightforwardly compute its long- 
term variation as the sum of two parts. The first one is 
purely Keplerian, and it vanishes over one orbital period Pb- 
The second one is due to the non-Keplerian variations of all 
the orbital elements induced by the dynamical perturbation 
considered. The total result is, thu^J, 



\ dt / 



AY 



C[ 



BY df djvl 
df dM dt 



+ 



(1) 



E«Hf] (§ )df 



In it, dM/dt and dn/dt are the instantaneous variation^ 2 ] 
of the Keplerian orbital elements computed with, e.g., the 
Gauss variation equations and evaluated onto the unper- 
turbed Keplerian ellipse, while df /dM and dt/df are the 
usual Keplerian expressions for such derivatives: see eq. {5J 
and eq. @ below. To make easier for the reader to go 
through the details of the following calculations, we list in 
Table [1] all the symbols used along with their definitions. 
More details are also given in the text. 

In order to give just some preliminary numerical es- 
timates of the effects computed, we will consider a rather 
typical star-planet scenario summarized in Table [21 It in- 

1 The analytical calculations have been performed with the aid 
of the software MATHEMATICA. 

2 Actually, dM/dt is the sum of the Keplerian mean motion n 
and a non-Keplerian term, as we will see later. Its Keplerian part 
yields from cq. {TJ the Keplerian variation of Y. 



Table 2. Reference stellar and planetary parameters adopted in 
the t ext. We use M* = Mp,, R* = R p,, J% = jf llPireaux et al.l 
120071) , S* = Sr ^Piiperslll99Sl , 120031 1. k^ = 0.03 jClarei| | l99^ 
m p = mj up , r p = rj up , fc2 p = 0.6 l lRagozzine fc W olf 2009). 
Concerning j%, we comp uted it from j% = (tep/3)(n 2 rp/Gm p ) 
llRagozzine fc Wolfll2009i ). where we assumed that the rotational 
frequency of the planet is equal to the orbital one n because 
of tidal effects. For simplicity, we assume that the equators of 
both the star and the planet coincide with the orbital plane, so 
that 'J'* = ipp = 0. Moreover, we assume an edge-on orbital 
configuration, i.e. i = 90 deg. 



Parameter (units) Numerical value 



M* (kg) 


1.98895 x 10 30 


R* (m) 


6.9599 x 10 8 


S* (kg m 2 s" 1 ) 


190.0 x 10 39 


Jl 


2 x 10~ 7 


£.'2* 


0.03 


m p (kg) 


1.89 x 10 27 


r p (m) 


6.991 x 10 7 


&2 P 


0.6 




3 x 10" 4 


a (au) 


0.04 


n (s" 1 ) 


2.5 x 10- 5 


Pb (d) 


2.9 


e 


0.07 


i (deg) 


90 


** (deg) 





V> P (deg) 






eludes a Sun-like host star harboring a close Jupiter-sized 
planet moving in the equatorial plane of the star along a 
moderately elliptic orbit. The orbital angular momentum 
is perpendicular to the line of sight, so that the orbit is 
in an edge-on configuration. We assume that the rotation 
of the planet has been locked by tidal effects to the or- 
bital frequency n, and that all the spins of the system are 
aligned. The meaning of the symbols used in Table [5] has 
been explained in Table [l] However, it must be noted that, 
if on the one hand, such a scenario is largely representa- 
tive of most of the transiting exoplanets discovered so far 
|Torres et al.ll2008l ). on the other hand some planets show- 
ing high inclinations of their orbits to the stellar equa- 
torial planes have been discovere d in recent times. They 
are WASP-33b dCollier Cameron et"ai] 120101) . WASP- 2b 



JCollier Cameron et al 



20071 ). WAs"p |5b J Anderson et all 



2008; iGillon et all |200£ ), W ASP-8b (iQueloz et al.l |2010|) 



WAS P- 15b dWest et al.ll2009l), WASP- 17b J Anderson et all 
120101) and HAT- P- 7b JNarita et al.ll2009l ; IWinn et al.ll2009l ) 



Our general expressions are valid also for such class of ob- 
jects. Moreover, the Kepler mission is expected to detect nu- 
merous planets with a wide distribution of both masses and 
semi-major axes; to a certain ex tent, it has already started 
to do that (|Borucki et al.ll2010l ). A quantitative extension 
of our results to such specific exoplanetary scenarios and to 
others exhibiting not yet explained features which may be 
accommodated by the effects investigated here may be the 
subject of a further, more applicative paper. 

Let us note that the validity of the formulas obtained 
is not limited only to exoplanetary scenarios. Indeed, they 
may be applied also for, e.g., designing new observational 
tests of general relativity in our solar system with natural 
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Table 1. Symbols used in the text along with their definitions. 

Symbol Definition 

G Newtonian constant of gravitation 

c Speed of light in vacuum 

M Mass of the primary in the two-body problem 

R e Equatorial radius of the primary in the two-body problem 

J2 Adimcnsional quadrupole mass moment of the primary in the two-body problem 

S Angular momentum of the primary in the two-body problem 

M* Mass of the star 

R* Radius of the star 

5* Angular momentum of the star 

J '2 Adimcnsional quadrupole mass moment of the star 

Ai2* Love number of the star 

m-p Mass of the planet 

r p Radius of the planet 

&2p Love number of the planet 

J2 Adimcnsional quadrupole mass moment of the planet 

lZ g = G(M+ + m p )/c 2 Gravitoelectric characteristic length of the system 

Xg = 2GS'*/c 2 Gravitomagnetic characteristic volume per unit time of the system 

Cx Tidal gravitational potential of a distant third body X 

rax Mass of the third body X 

rx Distance of the third body X from the star (assumed constant over one orbital period of the planet p) 

/Cx = Gmx/i 1 ! Tidal parameter of the third body 

lx = {lx, ly,lz} Unit vector pointing towards the third body (assumed constant over one orbital period of the planet p) 

r Star-planet relative distance 

a Semi-major axis of the orbit 

n = WG(Mi, + mp)/d A Keplerian mean motion of the orbit 

Pb = 2n/n Orbital period 

e Eccentricity of the orbit 

p = a(l — e 2 ) Semilatus rectum of the orbit 

v Speed of the planet along its orbit (assumed elliptical) 

i Inclination of the orbit to the plane of the sky 

7* Inclination of the angular momentum of the star to the plane of the sky 

'J'* Inclination of the orbit to the equatorial plane of the star 

ipp Inclination of the orbit to the equatorial plane of the planet 

/ True anomaly of the orbit 

u) Argument of pericenter of the orbit 

SI Longitude of the ascending node of the orbit 

u = u + / Argument of latitude of the orbit 

M Mean anomaly of the orbit 

{k} Set of Keplerian orbital elements (apart from the mean anomaly) 

Vg = Xg na ~ 2 (1 — e 2 )~ 7 ' 2 (l +ecos/) 3 Gravitomagnetic characteristic acceleration parameter of the system 

R Unit vector along the radial direction of the co-moving frame 

T Unit vector along the transverse direction of the co-moving frame 

N Unit vector along the normal direction of the co-moving frame 

A Perturbing acceleration caused by relativity, oblatcncss, etc. 

An = A ■ R Radial component of the perturbing acceleration 

At = A ■ T Transverse component of the perturbing acceleration 

Ajv = A ■ N Normal component of the perturbing acceleration 

S* Unit vector along the angular momentum of the star 

L Unit vector along the orbital angular momentum 

p Unit vector along the line of sight (pointing towards the observer) 

Y Generic observable quantity 

Y Instantaneous time variation of the generic observable quantity Y 

AY Net change of the generic observable quantity Y after one orbital revolution 

( Y ) Averaged time variation of the generic observable quantity Y 

t Observational time interval 

Uy Accuracy in measuring Y 

S (S) Latitude of the transit for elliptic orbits (circular orbits) 

Atj Transit duration 

V p Radial velocity 

K Semiamplitudc of the radial velocity curve 

Vo Radial velocity of the center of mass of the system 

At ec i Time span between the primary and the secondary eclipses 

Ptr Time span between two consecutive primary eclipses 
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and artificial bodies. In principle, they can be adapted to 
modified models of gravity as well. The latter aspect will be 
developed in further studies. 

The paper is organized as follows. In Section [2] we il- 
lustrate the main features of the classical and relativistic 
perturbing accelerations considered useful for our purposes. 
The time variations of the duration transit are worked out 
in Section [3] Section [4] deals with the time changes of the ra- 
dial velocity. Section[5]is devoted to the time interval elapsed 
between primary and secondary eclipses and its long-term 
temporal variations. The time span between consecutive pri- 
mary transits is investigated in Section [6] In Section [7] we 
summarize our results. 



Then, Ar, At, An are to be substituted into the right-hand- 
sides of the Gauss equations for the variations of the K eple- 
rian orbital elements. They are (|Rovll2005l ; ISoffelll 1989ft 



dt 



dt 



dll 

dt 




A R esmf + A T (z)], 

{ A« sin / + A T [cos / + ±(l -£)]}, 

=An (-) cosw, 



\ a / 



2 THE PERTURBING ACCELERATIONS 

Here we deal with a generic perturbing acceleration A in- 
duced by some classical and general relativistic dynamical 
effects. 

First, A has to be projected onto the radial, trans- 
verse and normal orthogonal unit vectors R, T, N of the co- 
moving frame of the test particle orbiting the central body 
acting as source of the gravita tional field. Their co mponents, 
in cartesian coordinates, are l|Montenbruck fc GillH2000l ) 



(cos fl cos u — cos * sin Q, sin u 
sin f2 cos u + cos * cos Q sin u 
sin * sin u 

(— sin u cos Q — cos * sin Q cos u 
— sin Q, sin u + cos * cos fl cos u 
sin * cos u 

(sin * sin Q \ 
- sin* cos Q 
cos* / 



(2) 



(3) 



(4) 



In eq. (J2}-eq. @, fi, u), * are the longitude of the ascending 
nodqj, the argument of pericenter, reckoned from the line 
of the nodeqj, and the inclination of the orbital plane to 
the reference {xy} plane, respectively. Moreover, u = f + 
uj is the argument of latitude. Subsequently, the projected 
components of A have to be evaluated onto the Keplerian 
ellipse 



1 + e cos / 



a(l 



(5) 



where p is the semilatus rectum and a, e are the semi- 
major axis and the eccentricity, respectively. The carte- 
sian coordinates of the K eplerian motion in space are 
l|Montenbruck fc GilllEfflOOft 



r (cos fl cos it — cos * sin Q sin u) 
r (sin SI cos u + cos * cos Q sin u) . 



(6) 



r sin* sin u. 



3 See also eq. l28t below. 

It is the intersection of the orbital plane with the equatorial 
plane. See eq. I I26II and eq. I I28II below. 



dui 

dt 



dM. 
dt 



= -cos*(f) + ^M«cos/+ 
+ At (l+ =) sin/] , 

_ Q-l± t_ A cog f+ 



+ At 1+S sin/ , 



(7) 

where n = \J ' G(M+ + m p )/a z is the Keplerian mean motion 
related to the orbital period by n = 2n/P^. 

As explained in the Introduction, the right-hand-sides 
of eq. (0, computed for the perturbing accelerations of the 
dynamical effect considered, have to be inserted into the 
analytic expression of the time variation dY/dt of the ob- 
servable Y of interest which, then, must be averaged over 
one orbital revolution according to eq. (JT| by means of (|Rovl 
12005ft 



df= (-Y Vl-e 2 dM, 



and 



dt 



(1 



=2\3/2 



n(l + ecos/) 



:df. 



(8) 



(9) 



2.1 The effect of the stellar oblateness 

The external gravitational field of a rotating body under- 
goes departures from spherical symmetry because of the dis- 
tortion of its shape due to the resulting centrifugal force. 
An oblate body of equatorial radius R e and adimensional 
quadrupole mass moment J2 affects the orbital motion of 
a test particle with a non-central perturbing acceleration 
ijCunn ingham 197dl : IVrbikll2005ft 



|(J2) 



3J 2 RjGM 

2r 4 



I [l - 5(r • 5) 2 1 r + 2(r • S)Sl\ 



(10) 

where S is the unit vector of the body's angular momentum, 
chosen here as z axis so that the equatorial plane is the ref- 
erence {xy}. We will adopt such a choice for the coordinate 
axes also for the other dynamical perturbations. According 
to eq. ©-eq. (@J and eq. (|6)l , the R — T — N components of 



cq. (fTO)l arc 



l(Jz) 



~8a(l-e 2 ) 4 
,2 



(1 + e cos/) 4 (1 + 3 cos 2*+ 



+ 6 sin 2 * cos 2m) 



-4 



A 



(J 2 ) 



"2a(l-e 2 ) 4 

3n 2 fl^J 2 
~2a(l-e 2 ) 4 



(1 + e cos /) 4 sin 2 * sin 2m, 



(1 + ecos/) sin 2^& sin -u. 



(11) 

Note that [r^R^a' 1 ] = L T~ 2 . For * = 0, i.e. for equato- 
rial orbits, only the radial component is not zero. For polar 
orbits, i.e. for "$> = 90 deg, the normal component vanishes, 
contrary to the radial and transverse ones. 



2.2 The effect of general relativity 

In its slow-motion and weak-field approximation, the Ein- 
stein's general theory of relativity predicts that a slowly 
rotating central body of mass M and proper angular mo- 
mentum S induces two kinds of small perturbations on the 
orbital motion of a test particle . The largest one is dubbed 
gravitoelectric ( Ma shhoonl 120071 ) . and depends only on the 
mass M of the body which acts as source of the gravitational 
field. It is responsible of the the well-known anomalous sec- 
ular precession of the perihelion of Mercury of 43.98 arcsec 
cty -1 in the field of the Sun. There is also a smaller pertur- 
bation, known as gravitomagnetic (|Mashhoonll2007l ). which 
depends on the angular mo mentum S of the central body: it 
causes the Lense-Thirring (L ense fc Thirrin e 1918) preces- 
sions of the node and pericenter of a test particle. 
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l|Soffellll989r i 

A<£ M) = %cos*(l + ecos/), 



.4 



(GM) 



-77 g ecos\l/sin/, 



^ M) = %sin*(l + ecos/)[2sinu+ 



(14) 



+ e 



sin / cos 
1 + e cos 



?)] 



with 



and 



^= a2( i-: T / 2 (i+ecos/)3 ' 



Xg 



2GS 



(15) 



(16) 



Note that [x 3 ] = L 3 T~\ so that [t] g ] - L T -2 . For equa- 



,(GM) 



(GM) 



.(GM) 



torial orbits ,4^ lvl; = and A^ 1 ' / 0, A£* m; / 0. Instead, 
for polar orbits, i.e. for ^ = 90 deg, only the normal com- 
ponent does not vanish. 



2.3 The tidal bulges 

By neglecting the lag due to dissipation w hich is negligible 
for giant planets (jMurrav fc Dermottll2000r ). the perturbing 
acceleration due to the tidal bulge rais ed on the plan et by 
the host star is entirely radial, so that l|Cowling||l938h 



.4 



(tid P ) 



3fc 2p r|G(M,+m p ) 



A 



(tid P ) 



(17) 



2.2.1 The gravitoelectric, Schwarzschild-like perturbation 
By defining 

G(M* + m p ) 



Ha 



(12) 



where c is the speed of light in vacuum, the R — T—N com- 
ponents of the ge neral relativ istic gravitoelectric perturbing 
acceleration are (|Soffell ll989) 



A 



(GE) 



l^fy (1 + e cos /) 2 (3 + 2e cos / - e 2 + 



+ 4e 2 sin 2 /) 



,(GE) 



A^' = 7^M + ecos/) 2 4esin/(l + e cos/), 



.4 



(tid P ) 

AT 



0. 



Here A?2 P is the Love number of the planelQ It measures 
how the mass redistribution induced by the non-uniform 
stellar potential actually affects the external gravity field 
of the planet. Its values are in the range kz p ss 0.1 — 0.6 
(|Ragozzine fc Woldl2009l ). 

A similar expression holds for the effect due to the tidal 
bulge raised by the planet on the star: the required substitu- 
tions are M* — > m p ,m p — > M*,r v — > 7?*, fc2 P — > fe*. Bodies 
with most of their mass near their cores, like main-sequence 
stars, have very low fe*. It is so because the distorted outer 
envelope has little mass, t hus affecting negligibly the outer 
gravitational field. Indeed. IClaretl (|1995T ) yields &2* ~ 0.03. 



,(GE) 



0. 



Note that [Tl g ] = L, so that [n 2 TZ g ] = LT 



(13) 



2.2.2 The gravitomagnetic, Lense-Thirring-like 
perturbation 

The R — T — N components of the general relativistic grav- 
itomagnetic perturbing acceleration induced by the rotation 
of the central body with proper angular momentum S are 



2.4 A distant third body X 

The perturbing potent ial induced by a distant, third body 



the perturbing potent i. 
Xis l|Hogg et al.lll99ll ) 



*-^[r>-3Mx)'] 



24 



(18) 



5 Here &2 P is t wice the Love n umber customarily used in binary 
stars literature (Cowling 1938). 
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where lx = rx/rx is a unit vector pointing towards X. By 
denoting l x ,l y ,l z the direction cosines of rx, i.e. the com- 
ponents of lx, we will express Ux as 

Ux = ^- [r 2 - 3(xl x + yl y + zl z f] . (19) 

After working out the cartesian components of the perturb- 
ing acceleration Ax = — VfJx from eq. (|19[) . the explicit 
expressions of the R — T — N components can be worked out 
with the aid of eq. ©-eq. Q and eq. ([6]). They are 

4?' = f rfe^[*.(**,n,fc)-i»u+ 



+ 91 



( X ) 3a(l- e 2 )>Cx 



I 2 
:(£7,Z X )cOSU , 



(1 + e cos /) 



£7, lx) sin2w+ 



£ c (* *,€l,lx) cos 2ul 



(20) 



A 



where 



(X) 3a(l- e 2 )?Cx 

(1 + e cos /) 



[«n s (**,n,fx)sin«+ 



+ *Jt c (**,£M x )cosuj , 
Gmx 



k\ 



(21) 



is the so-called tidal parameter; [fCx] — T~ 2 . Note that both 
lx, which is assumed constant over one orbital revolution of 
the perturbed planet, and £2 here refer to the stellar equa- 
torial plane. 

The coefficients of eq. (|20[) are 



$R S = l z sin** + cos*]/* (l y cos £7 — l x sin £7) , 
£Hc = lx cos £7 + l y sin £7, 



Is = 



Z s sin\E'*+cos\E'*(ZyCosn — Z^sinn) — (Z^cosfi + Zysinfil 



Z c = (^ cos S1 + ly sin £7) [J z sin ** + cos ** (l y cos £7— 

— J^sinfl)], 

Vis = [l z sin** + cos** (l y cos £2 — l x sin £7)] x 

x [l z cos ** + sin ** (Ja: sin £1 — l y cos £7)] , 

9tc = (Ja; cos £7 + Z y sin £7) [Z z cos ** + sin** (l x sin £7— 

— l y cos £2)] . 

(22) 
Note that, according to eq. (|22[) . eq. (|20[) vanishes neither for 
** = nor for l z — 0. If both the planet and the perturber 
X are coplanar and lie in the equatorial plane of the star, 
then A™ = 0, while Af ] / 0, Af ] / 0,Af } / A^\ as 
intuitively expected. 

Actually, in view of their complexity, the previous for- 
mulas do not allow one to obtain manageable analytic ex- 
pressions of general validity for the averaged time variations 
of the observable quantities we are interested in. Thus, in the 



following, we will release approximate expressions for such 
kind of perturbation. However, the exact effects of X can 
be numerically computed in specific arbitrary exoplanetary 
scenarios following the procedure discussed in Section [T] 



3 TIME VARIATIONS OF THE TRANSIT 
DURATION IN ELLIPTIC ORBITS 

For an unperturbed Keplerian elliptic orbit arbitrarily in- 
clined to t he line of sight, the tran s it duration Atd can b e 
written as (|Tinglev fc Sacketdl200 j ; Ijordan fc Bakos]|2008l ) 



A+ 2(7?* + r p ) 

Atd — — — cos o. 



(23) 



In it, Ri, and r p are the radii of the star and the pl anet, 
respe ctively. The Keplerian speed of the planet v is (|Rovl 
120051 ) 



VT 



:\A + 2ecos/ + , 



(24) 



The latitude S of the transit on the stellar disk is defined 
from 



sin 8 = 



R* + r p ' 



(25) 



in which r is given by eq. (|5|. The parameter i is the angle 
between the unit vector L of the orbital angular momentum 
and the unit vector p of the line of sight poi nting towards 
the observer. From the spherical law of cosines (JGellert et al.l 
ll989l ; IZwillingerlll995l ) 

cos B = sin C sin A cos b — cos C cos A (26) 

with the identification^] A — >• **, 73 — >• tv — i.C — >• I*,b — ► 
7r — £7, it turns out 



cos i = sin ** sin 7* cos £2 + cos ** cos 7* 



(27) 



where ** is the angle between L and the unit vector S* of 
the star's proper angular momentum, 7* is the angle between 
»§* and p, and £1 is the longitude of the ascending node 
defined from 



sin ** sin 7* cos £7 — (S* X L) ■ (S 1 * X p). 



(28) 



Incidentally, it may be noted from eq. (|27| that for equato- 
rial orbits, i.e. for ** = 0, we have i — 7*. Thus, Atd, for an 
elliptic orbit, is a function ofj f,a,e through v and r in cos 5, 
and of ** and £7 through i in cos S. As a consequence, when 
the effect of a given dynamical perturbation on Atd is con- 
sidered, the variations of all such Keplerian orbital elements 
have to be fully taken into account. 

In the following calculations of the time variation of Atd 
due to various dynamical perturbations the expression 



a(l 



(7?*+r p )(l + ecos/) 



-1/2 



(29) 



6 In such a way, Q results to be prograde with respect to the 
orbital motion, i.e. £2 follows it, coherently with the definition of 
the longitude of ascending node. 

7 We neglect possible time changes of the star's spin axes. 
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appears: to make computation simpler, we will approximate 
it with 



a(l — e 2 ) cosi 



(R* + r P ) 



-1/2 



(30) 



3.1 The effect of the stellar oblateness 

An exact calculation, valid to all orders in e, cannot be 
performed because of the difficulty of the integrals involved 
in the evaluation of eq. {l}. Thus, in addition to eq. (|30[l . we 
have to approximate also eq. (|24|l entering the integrands 
with negative powers. We will retain the first order terms in 
their expansions, which is admissible for e < y/2 — 1 ss 0.4. 
Such an approximate calculation shows that the over- 
all averaged time variation of At d induced by the stellar 
oblateness does not vanish for elliptic orbits. It turns out 
that the non-vanishing terms are those due to the pertur- 
bations in o, e, \t*, Cl, u),M. This is a clear example of how 
simply inserting the long-term precessions of the Keplerian 
orbital elements into instantaneous expressions would yield 
incorrect results; indeed, the long-term variations of a, e, *&* 
caused by the oblateness of the primary are notoriously zero 
JRovl [2005V By defining 



sin 8 



(R* + r P ) 



(31) 



we finally have 
|(JJ) 



dt / \ a I 



K,\ 2 3J|i 



;in 2M* * sin 7* sinfi 



2 v / l-(l-e 2 ) 2 sin 2 S 



- + 



+ 0(e 2 ). 
In the limit e — > 0, eq. (|32|1 reduces to 
dAt d 



dt 



(J 2 *) 



2 tan 8 sin *)/* sin 7* sin D. / dil 

dt 



(J 2 *) 



(33) 
obtai ned in lloriol (|2010l ) for circular orbits. Indeed (JRovl 

Hool), 



( J 2*)\ Q ( T> \2 

= -\ nJ i\^) cos ** (34) 



~~dt 



for e — > 0. For another derivation of eq. (|33p . which can be 
applied to all perturbations when e — > 0, see Section 13.41 
Note that, for exactly edge-on orbits, both eq. (|32p and eq. 
(|33[) vanish. The same occurs if the orbit is equatorial, ir- 
respectively of its inclination to the plane of the sky. It is 
important to note that the long-term effects of eq. (|32[) -eq. 
(|33[) are not secular trends linearly changing in time because 
of the presence of ^* and Q as arguments of trigonomet- 
ric functions. Indeed, such Keplerian orbital elements do, in 
general, slowly vary in time because of dynamical effects like 
a distant third body X, general relativity, and the oblateness 
itself. 

Finally, let us mention that the oblateness jf of the 
planet also affects Atd in the same way. The related formulas 
can be obtained from the previous one provided that the 
substitution _R* - 



J2 — > J2 is performed. 



3.2 The effect of general relativity 

3.2.1 The gravitoelectric, Schwarzschild-like perturbation 

From eq. (|13|) it turns out that the general relativistic gravi- 
toelectric term can only affect dAtd/dt through its perturba- 
tions on a, e, M. The long -term variat ions for such Keplerian 
orbital el ements are zero ((Soffel|[l989J) , apart from the mean 
anomaly (|loridl2007J V 

Instead, after an approximate calculation valid for e < 
^/2 — 1 « 0.4 analogous to the previous one, it turns out 
that 



dAt d 



dt 



<OE) 



0. 



(35) 



The same holds also for circular orbits, as can be im- 
mediately seen by inspecting the general method of Section 
I3.4l which involves the long-term variations of ^* and Q. 

3.2.2 The gravitomagnetic, Lense-Thirring-like 
perturbation 

According to an approximate calculation in e, valid for 
e < \/2 — 1 ~ 0.4, analogous to the previous ones, the aver- 
aged time variation of Atd induced by the general relativistic 
gravitomagnetic effect is, in general, non-zero for elliptic or- 
bits. In particular, the non- vanishing terms are those due to 
the perturbations in ^* and £1. This is another example of 
how the mere insertion of the long-term precessions of the 
Keplerian orbital elements into instantaneous expressions 
would yield incorrect results; inde ed, the Lense-Thirring 
long-term p recession of <&* vanishes (|Lense fc Thirringlll918l ; 
ISoffe]||l989r ). 
We have 



(32) / dAt d 



W GM) \ 

dt \ I 



- < " S * ) 

\ c' 2 a' A n l 



sin S sin gf * 



2 V / l-(l-<= 2 ) 2 



x [£(e,n,w,/*,**) + 9Jt(e,n,w J /* 



(36) 



with 



£ = e 2 (1-e 4 ) (cos^sin**- 

— cos ^1/* sin/* cos fl) sin 2ui, 

M = {I - e 2 ) U\l - e 2 (l - \ 

— e (l + e ) cos2o;} sin/* sinSl; 



(37) 



the term containing £ comes from the perturbation in ^*, 

while the term with *Xt is due to dQ/dt. 

For e -> 0, eq. (gHJ-eq. §7} yield (|lorioll2010h 



dAt, 

dt 



|(GM) 

(I 



2 tan 5 sin SI* ^ 



1/* sinQ (2GS 



2 tan i5 sin ty ^ sin I \ i 



/2GS* \ _ 
\S*aT) - 



d£2|(GM)\ 
dt I / ' 

(38) 
According to eq. (|36[) -eq. (|38p . also the long-term grav- 
itomagnetic time variation of Atd is a harmonic one. For 
exactly edge-on orbits, both eq. Q36p an d eq. (|38p vanish. 
The same occurs if the orbit is equatorial, independent of 
its orientation with respect to the line of sight. 



8 L. Iorio 



3.3 The tidal bulges 

It is found that the long-term variations of the duration 
transit caused by the tidal bulges mutually raised by the 
planet and the star on each other vanish. 



3.4 A distant third body X 

It turns out that a distant body X induces a long-term har- 
monic variation of the duration transit which does not vanish 
in the limit e — > 0. In this case, we have 



dAt a 



dt 



(x) 



^)3tan53E(fx,**,i*,n), (39) 



a distant body X, it turns out from eq. (0 
= 0, 



dal(X) 
dt I 



dM>* |(X)\ 

dt I / 



dOj(X) 
dt I 



= 3Kx( 2 1 ~ e2)3 (lx cos tt + l y sin fi) x 

x [l z cos ** + sin ** (l x sin Q, — l y cos SI)] , 

= 2^ 1^ Sm ** "*" COS ** X 

x (ly cos SI — l x sin SI)] [l z cos**+ 

+ sin ** (l x sin Q, — l y cos SI)] , 

(44) 
contrary to the stellar oblateness and the general relativistic 
gravitomagnetism for which the long-term precessions of ** 
vanish. It turns out that eq. (|44[) . substituted in eq. (|42[) -eq. 
((43)) . yields the same result of eq. j39])-eq. (|40|) . 



with 

X = [l z cos*l/* + sin** (— l v cos SI + l x sinfi)] x 

x {-/z cos** sin/* + sin** [l z sin/*sinn+ (40) 

+ cos/* (l x cosO + ly sin SI)]} . 

Note that if both the perturber body X and the perturbed 
planet lie in the equatorial plane of the host star, i.e. if 
** = and l z = 0, then 



dAt„ 



dt 



(X) 



= 0, 



(41) 



independently of the inclination of the orbital plane to the 
plane of the sky. Instead, if ** = 0, but X is not coplanar 
with the perturbed planet, the induced long-term variation 
of the duration transit does not vanish. The duration transit 
is not affected by X if the orbit of the perturbed planet is 
exactly edge-on. 

It is interesting to note that eq. (|39p . with eq. (|40p . 
can be also obtained in another, simpler way; the following 
approach is valid for all the dynamical perturbations pre- 
viously examined. Indeed, for circular orbits the expression 
for Atd of eq. (|23f) gets simplified in such a way that 



3.5 Numerical evaluations 

A feature common to all the variations of Atd previously 
examined is that they all vanish for exactly edge-on orbits. 
Moreover, independently of the inclination of the orbital 
plane to the line of sight, they vanish also for equatorial 
orbits, apart from the case of X unless it is coplanar with 
the perturbed planet moving in the star's equatorial plane. 
Thus, we assume in the following i = 87 deg and ** = 
15 deg. With such assumptions we have 



dt 



dAf,, 
dt 



dAf„ 
dt 



dAf,, 
dt 



dAf„ 
dt 



dAf,, 
dt 



<j; i 



(XI) 



(GM) 



(tid *) 



(tid p) 



(GE) 



< 1 x 10" 



< 9 x 10~ 10 , 



< 1 x 10" 



= 



(45) 



= 0. 



dAt„ 



dt 



2 tan S sin i I di\ dAtd (da 
n \dt J da \ dt 



(42) 



From eq. (|27|) it turns out 



sini (5t) = (sin ** cos 7* - cos** sin/* cos 17) ^§f + 

+ sin** sin/* sin SI (^-) . 

(43) 
Thus, only for circular orbits it is admissible to straightfor- 
wardly use the long-term variations of a, ** and SI, as shown 
by eq. (l42[) -eq. (|43[) . In the case of the perturbation due to 



By assuming e = 0, a third body X with the mass of Jupiter 
and located at rx = 0.1 — 1 au yields 



dAt a 



dt 



(Jup) 



<8x 10" 



8 x 10" 



(46) 



according to 1391 The effect of an Earth-like perturber rx 
0.1 — 1 au is, instead, 



dAt c 
dt 



(Ear) 



< 2 x 10~ 7 - 2 x 10" 



(47) 



The time variations of the duration transit induced by 
the planetary and stellar oblateness for closer orbits, i.e. for 
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Duration transit: stellar quadrupole mass moment 
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Figure 1. Maximum values of the long-term time variations ( At d \ as a function of a (0.014 au < a < 0.449 au) for different values of 
the eccentricity: e = 0.005 (red dash-dotted line), e = 0.03 (blue dashed line), e = 0.1 (light blue dotted line), e = 0.4 (yellow continuous 
line). For the stellar and planetary physical parameters we used the standard values of Table [2] We adopted i = 87 deg, <P* = 15 deg, 
while we fixed the periastron at u) = 90 deg; we also considered sin/* sinQ m 1. 
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Figure 2. Maximum values of the long-term time variations ( At<j > as a function of a (0.014 au < a < 0.449 au) for different values of 
the eccentricity: e = 0.005 (red dash-dotted line), e = 0.03 (blue dashed line), e = 0.1 (light blue dotted line), e = 0.4 (yellow continuous 
line). For the stellar and planetary physical parameters we used the standard values of Table [2] We adopted i = 87 deg, >t* = 89 deg, 
while we fixed the periastron at ui = 90 deg; we also considered sin/* sinQ rs 1. 



a — 0.015 au, would be of the order of 7 x 10 7 (j|) and 2 x 
10 -9 ( J| ), respectively. In Figure[l]we plot the upper bounds 

in (Atd) due to all the dynamical effects considered, apart 
from the third body X, as a function of a for different values 
of the eccentricity. We let a vary within the range of the 
so-far discovered transiting exoplanetqj, with a maximum 
value of e yielding periastron distances not smaller than the 
star's radius. We adopted nearly edge-on (i — 87 deg) and 
almost equatorial (vp* = 15 deg) orbital configurations. The 
physical parameters of the star and the planet have been 
retrieved from Tabled It is interesting to look also at almost 
polar orbital configurations ('J* « 90 deg). This is done in 
Figure[2]obtained for the same values of Figure[l] apart from 
the inclination of the orbital plane to the stellar equator set 
now to <]/* = 89 deg. 



3.6 The observability of the time variation of the 
transit duration 

The accuracy obtainable in measuring the transit duration 
is of the o r der of 1.5 — 5 s; see the accurate discussion in 
iFord et all <|2008r i and IJordan fc Bakosl (|2008l ). Thus, over 
an observational time span of r = 10 yr it should be possible 
to achieve an accuracy in measuring the time variation of the 
transit duration of the order of 



5 x 10" 



1.6 x 10" 



(48) 



A similar conclusion is reached by iMiralda-Escudel ((.2002). 
Thus, it would be very difficult to measure the effects pre- 
viously examined with the time variations of the duration 
transit in typical exoplanetary scenarios, with the possible 
exception of a third body X. 



4 THE RADIAL VELOCITY 

The basic observable in spectroscopic studies of exoplanets, 
transiting or not, is the radial velocity V p . Its expression for 
elliptic orbi ts, up to the velocity of the system's center of 
mass Vo, is <|Battenll2o6~th 



V p — K [e cos uj + cos(/ + w)] . 



(49) 



where 2K is the total observed range of radial velocity de- 
fined by 



K 



VT 



(50) 



See on the WEB |http://www,exoplanet.eu| 



Perturbing dynamical effects affect the radial velocity 
as well by inducing, in principle, a non-vanishing net radial 
acceleration over one orbital period. It can straightforwardly 
be worked out from eq. JTJ) with Y — > V p by noting that, in 
this case, the perturbations of all the six Keplerian orbital 
elements are involved. 



4.1 The effect of the stellar oblateness 

The stellar oblateness causes a long-term harmonic variation 
of the radial velocity only if the orbit is elliptic. Indeed, it 
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turns out 



K 



u 2 *) 



32(l-e 2 ) 7/2 sini' 

x [J"c(e, i,i*,ft, 'J'*) cosaj+ 
+ Js(e, i, I*,u, ft, **)sinw] 



(51) 



with 



J c = 10(1 - e 2 ) cos i sin/* sin 2** sin ft, 

J" s = 2(1- e 2 ) cos j sin 2** (cos /* sin **- 

- sin /* cos ** cos ft) + sin 2 i [7 + 47 cos 2*+ 
+ sin 2 ** cos 2lj - f| (259 + 429 cos 2**- 

- 44 sin 2 **cos2a;)] . 

(52) 
It is an exact result in e, and vanishes in the limit e — > 0. 
Note that, according to eq. (|52p . eq. ()51[) vanishes neither 
for equatorial orbits nor for edge-on configurations. 



4.2 The effect of general relativity 

Also general relativity affects the radial velocity of non- 
circular orbits. Indeed, an exact calculation in e yields the 
following long-term harmonic signatures. 



KGE)\ 



r=(»%) 



15e(l + e )sinisino; 
8(l-e 2 ) 5/2 ' 



(53) 



and 



V t 



(GM)\ 



, / = (^) irT ^ F [V c (i,/*,fi,**)cos^+ 



(54) 



+ V s (i,/*,ft, **)sinw], 

with 
V c = 11 cot i sin 7* sin** sin £2, 

V s = 2fi {cos ft sin 27* (sin** -sin 3**)- 

- sin 2 ** cos ** [cos 21* (3 + cos 2£2) + 2 sin 2 ft] + 

+ 104 sin 2 i cos**} . 

(55) 
For e — >• both eq. ()53[) and eq. (|54l) vanish. Also the general 
relativistic effects are non-vanishing either for edge-on or 
equatorial orbits, with eq. (|53|) which is independent of "I/*, 
contrary to eq. (154|l . 



4.3 The tidal bulges 

The tidal bulge of the planet induces a non-zero long-term 
harmonic variation of the radial velocity. Its exact expression 



(tid p) 



Vf 



3e[-496 + 5e 2 (-80 + 81e 2 +8e 4 )]fc 2 p sin i sin u 
X 128(l-e 2 ) 13 / 2 ■ 

(56) 
The long-term variation of the radial velocity due to the 
tidal bulge raised by the planet on the star is 



(tid *) 



vr 



3e[-496+5e 2 (-80+81e 2 +8e 4 )]fc 2 * sin i sin u 
X 128(l-e 2 )13/ 2 ' 

(57) 
They vanish for e — > and are independent of **. 

It is easy to understand that, concerning the tidal dis- 
tortions, those experienced by the planet are more effective 
than those suffered by its hosting star. Indeed, eq. (|56[) -eq. 
((57)) tell us that their ratio goes as 






a. 



fc2p 
fe*' 



(58) 



According to Table [21 the ratio of the star-planet masses is 
typically of the order of 10 3 , while the planet-to-star ratios 
for the radii and the Love numbers are about 0.1 and 20, 
respectively. 



4.4 A distant third body X 

A net non-zero long-term harmonic effect on V p is induced by 
a distant third body X. Its general expression for arbitrary 
values of the eccentricity and of the position of X is too 
cumbersome to be explicitly displayed here; it turns out 



vr) oc 



•wi-f> 2 W 2 
■(oACx) l Q . / gP(fl,a;,tt»,lx). (59) 



The quite complicated function 2U does not vanish for {„ = 
0, ** = 0, i.e. if both the perturbed and the perturbing 
planets lie in the star's equatorial plane. As for the stellar 
oblateness and general relativity, (Vp } vanishes for e — > 0. 



/ 
It should not be considered as contradictory because it is 

the otherwise constant velocity Vo of the system's center of 

mass that is altered by the perturber X, independently of 

the relative motion of the perturbed planet and the star. 
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4.5 Numerical evaluations 

The "standard" orbital scenario of Table [2] yields 

V p (tidp) )\ J x II! 7 tu , -'. 



|/t/ p (ge) \| < lxr'ms" 



Vp 



0?) 



v, 



(J$) 



(tid *) 



vf 



v P (GM) )\ 



< 6 x 10~ s m s~" 2 , 

< 3 x 10" 9 m s" 2 , 

< 2 x 10~ 9 m s -2 , 

< 9 x KT 11 ms -2 



(60) 



Concerning a distant third body X, for rax = mj up ,rx = 
0.1 — 1 au and l z — 0, i.e by assuming coplanarity with the 
perturbed planet, we have 



v; 



(.Tup) 



< 1 x 10 _b - 1 x 10" 



For an Earth-sized perturber at 0.1 — 1 au we have 



K (Ear) ^ <5x 10- 8 -5xl0" 



(61) 



(62) 

In Figure [3] we plot the magnitude of the upper bound 
of ( V p ) for all the dynamical effects considered, apart from 
a third body X, as a function of the semi-major axis for 
different values of the eccentricity. Also in this case, we let 
a vary within the range of the so-far discovered transiting 
exoplanet^j, with a maximum value of e yielding periastron 
distances not smaller than the star's radius. Table [2] has 
been used for the values of the physical parameters of the 
star and the planet, and for the various orbital inclinations 
i, **, ip p ; the periastron has been set equal to u) = 90 deg. 



4.6 Measurability of the long-term radial 
acceleration 

The present-day level of accuracy in measuring the radial 
velocity is about 2 m s , so that over r = 10 yr a net 
radial acceleration should be detectable at a ~ 3 x 10 -9 m 
s -2 level. Such a guess is confirmed by the case ofjfj HD 
126614A, monitored for the past 10 yr, for which a secular 
change of the radial velocity ([Howard et al.ll2010l ) 



V, 



(mcas) 



(5.13 ±0.06) x 10" 



(63) 



9 See on the WEB http://www.exoplanet.eu 

10 In addition to the planet HD 126614Ab (a = 2.35 au, e = 0.41, 
P b = 3.41 yr), the star HD 126614A has also a faint M dwarf com- 
panion of M(j w ?s 0.3Mq at a (sky-projected) distance of 33 au. 
Named HD 126614B, it was discovered with direct observations 
using ada ptive optics and the PHARO camera at Palomar Ob- 
servatory l|Howard et al.ll2010f ). An even more distant M dwarf, 
dubbed NLTT37349 and with Af dw sa 0.2M©, is also present 
in the system at a (sky-p rojected) distance of about 3 x 10 3 au 
iGould fc Chanamejliooj) . 



has been measured. Incidentally, let us note that for HD 
126614A none of the "internal" star-planet dynamical effects 
considered so far (tidal bulges, oblateness, general relativity) 
can explain 11 ! eq. ()63[) ; indeed, the largest one, caused by 
the general relativistic gravitoelectric field, amounts to just 
10~ 12 m s~ 2 . Concerning a distant third body X, it turns 
out that a Jupiter-sized planet at 2.2 au (Pb = 3.0 yr), 
assumed coplanar with HD 126614Ab in the equatorial plane 
of their hosting star, would cause a radial acceleration of the 
order of 5 x 10~ 7 m s~ 2 . Alternatively, a rocky Earth-like 
body at 0.32 au (Pb = 0.17 yr) or a Neptune-type planet 
(mx = 14m®) at 0.75 au (Pb = 0.6 yr) would induce the 
same effect. It must be recalled that the quality of t he fit 
of the 10 yr-long data record bv lHoward et al.l i|2010l ) was 
not improved, at a statistically significant level, by a two- 
planet model for a variety of masses and distances for the 
second one. A brown dwarf-like object (mx = 80raj up = 
7.6 x 10" 2 M Q ) at 9 au (P b = 24.4 yr) would also be a viable 
candidate. According to eq. flSSJl, neither NLTT37349 nor 
HD 126614B can accommodate eq. (|63[1 . Indeed, the effect 
due to the closest M dwarf would be as large as < 10 -8 m 
s~ 2 , while that due to the farthest one is as small as 10~ 14 

„, „-2 



5 THE TIME VARIATION OF THE 
TEMPORAL INTERVAL BETWEEN 
PRIMARY AND SECONDARY TRANSIT 

Another directly observable quantity in transiting exoplan- 
ets is the time elapsed Ai cc i between the primary and the 
secondary eclipses. The primary eclipse occurs when the ex- 
oplanet starts to obscure the host star, while the beginning 
of the secondary eclipse is when the planet starts to be oc- 
culted by the star. For circular orbits, At ec \ is simply half the 
orbital period, while for eccentric orbits we have a mo re com- 
plicated expression involving e and uj as well. It is (|Sterng 
1 19401 ; I Jordan fc Bakosll2008l ) 



Ai cc i — (taecl — tied) Y — ~^- 



- e cos id 



- (e sin uj) 



■ + 



+ arctan 



up to terms proportional to cot i. 



(64) 



5.1 Secular time variations of Ai cc i 

Non-Keplerian orbital perturbations affect, in principle, 
Ai cc i by making it vary orbit after orbit with a characteris- 
tic long-term pattern. It is possible to analytically work out 
exact expressions for all the dynamical effects considered, 
apart from the case of the distant third body X for which 
an approximate formula in e is released. Such long-term har- 
monic effects are 



|W) 



_R t \2 3eJ2 (3+5 cos 2*,) sinuj 



(¥) 



(65) 



11 We assume i = 90 deg. 
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Figure 3. Maximum values of the long-term time variation (Valiums 2 , as a function of a (0.014 au < a < 0.449 au) for different 
values of the eccentricity: e = 0.005 (red dash-dotted line), e = 0.03 (blue dashed line), e = 0.1 (light blue dotted line), e = 0.6 (yellow 
continuous line). For the stellar and planetary physical parameters, and for the inclination of the orbit to the plane of the sky and to the 
star/planet equators we used the standard values of Table [2] We fixed the periastron at lo = 90 deg. 



|U?) 



(GE) 



-r p \ 2 3ejJ(3+5cos2i/i p ) sinu . . . 

~ a > 2V'l-e 2 (l-e 2 sm2„) 2 ' ^ ' 



TZ a \ 12e^/l — e 2 sinu 



(67) 



l(tid P ) 



\ m p J V a J 



15e(-128+1152e 2 +6048e 4 +4200e 6 +441e s )k 2 psmaj 
64(l-e 2 )V2(l_e2 sin 2 [J ) 2 ' 

(69) 



/ dAt ccl K GM )\ _ / GS„ \ 24ecos^^ sin u 

V dt J \t?na*) (l_ e 2 s in 2 u) 2 



(68) 
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dAt ccl 
dt 



|(tid + ) 



(X) 



(l- e 2) 3 / 2 (l- e 2 ain 2„)3 

ES=i(W7)*®k (a «,**)]+ (71) 



+ 0(e 2 ), 

where (hlj) k , k — 1 ... 6 are to be intended as the six prod- 
ucts lx,lxly,lxlz,ly,l y lz,l 2 of the components of ix, while 
each of the six &k,k = 1...6 is a linear combination of 
trigonometric functions whose arguments are, in turn, linear 
combinations of Q, w, ^*. Note that eq. (|65|)-eq. (|7ip vanish 
in the limit e — > 0, as expected, while they do not vanish for 
** =0. 

It is interesting to note that the same considerations of 
Section 14.31 concerning the relative strengths of the tidally- 
induced effects on the radial velocity hold also in this case. 
Indeed, according to eq. H69[) -eq. (|70p the ratio of the planet- 
to-star tidal effects is equal just to eq. (|58H • 



5.2 Numerical evaluations 

For the typical star- planet scenario of Table [2] we have 

3 x 10" 7 , 



/ dAt ocl 
\ dt 


(tid p) 


> 


1 dAi ocl 
\ dt 


(GE)\ 




1 dAt ocl 
\ dt 


C£)\ 




1 dAt ocl 
\ dt 


Ca)\ 




1 dAt ocl 
\ dt 


(tid *)\ 


1 dAt ocl 
\ dt 


(GM)\ 





< 2 x 10" 7 , 



< 3 x 10" 



< 2 x 10" 



< 1 x 10" 



< 4 x 10" 



(72) 



The effect of a distant perturbing body X is comparable to, 
or even larger than eq. ([72l . Indeed, for rax = mj up ,rx = 
0.1 — 1 au and l z = 0, i.e by assuming coplanarity with the 
perturbed planet, we have 



dAt r 



dt 



(Jup) 



< 1 x 10" 



1 x 10" 



(73) 



If we assume an Earth-sized perturbing body X at 0.1 — 1 
au we have 



dAt c 



dt 



(Ear) 



< 4 x 10" 



4 x 10" 



(74) 



'"p \ (E±V 

M, ) \ a I 

15e(-128+1152e 2 +6048e 4 +4200e 6 +441e s )fc 2 * 
64(l-e 2 ) 7 / 2 (l- e 2 sin 2 ^) 2 

(70) 



Figure [4] illustrates the magnitude of the upper bound 
of ( Ateci) for all the dynamical effects considered, with the 
exception of X, as a function of the semi-major axis for dif- 
ferent values of the eccentricity. 



~5'.3 The measurability accuracy 

Conce rning the measurability of At e cl> IJordan fc Bakosl 
(2008) point out that it would be dominated by the un- 
certainties in the time of the secondary eclipse feed which is 

;80 s. 



evaluated by them to be of the order of about at 2 
Thus, it is reasonable to put 



0"At„ 



(J I., 



% tr 3 / 2 : 



(75) 



where 9l2tr is the number of secondary transits observed. By 
posing t = 10 yr and assuming VUiti = 10 3 it is possible to 
obtain 

^A tccl ^8xl0- 12 . (76) 

Thus, all the effects of eq. (172[) would be measurable. 



6 THE VARIATION OF THE TRANSIT 
PERIOD IN ELLIPTIC ORBITS 

Another important directly measurable quantity in tran- 
siting exoplanets is the temporal interval P tr be- 



tween successive primary transits (|Mi ralda-Escude 2002; 
Agol et all 120051; iNesvornv fc Morbidellil 120081: iNesvornvl 



20091; IVeras fc Fordl 120091 ; iPavne et alJl201Ch . To compute 
the variation of the transit period P tr in elliptic orbits we 
will make use of the so-ca lled guiding center approximation 
IjMurrav fc Dermottll2000h . In it, the Keplerian motion of a 
test particle p moving along an elliptic orbit with focus F 
is described in a suitabl e rotating reference frame {x y }. 
Referring to Figure 2.8 of lMurrav fc Dermottl (|2000h for the 
orientation of the axes, the guiding center G is the origin 
of such a frame. It moves with angular speed equal to the 
mean motion n along a circle centered in F ad having radius 
equal to the semi-major axis a. By orienting the y axis tan- 
gentially the guiding center's circle towards the motion of 
the test particle and the x axis perpendicula rly to it from 
F to p, it turns out (JMurrav fc Dermottl I2OO0T ) 



= rcos(f-M) 
= rsin(/-.M). 



(77) 



Concerning the argument of the harmonic functions in eq. 
(|77[1 . the exact rel ation between f and M can be derived 
from eq. ©. It is i|Capderoull2"ooll 



M = 2arctan 



1-e, // 

tan - 

1 + e V2 



e\J\ — e 2 sin / 
(1 + ecos/) 

A useful approximation of eq. (|78p is (Roy 2005) 
M 



f - 2e sin / + -e 2 sin 2/ + C(e 3 ) 



(78) 



(79) 



As far as the transit period P tr is concerned, from the choice 
of the axes in the rotating frame it turns out that the com- 
ponent relevant to such a phenomenon is y . Thus, the de- 
parture of Ptr from the Keplerian orbital period Pb = 2n/n 
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Figure 4. Maximum values of the long-term time variations ( At cc i \ as a function of a (0.014 au < a < 0.449 au) for different values of 
the eccentricity: e = 0.005 (red dash-dotted line), e = 0.03 (blue dashed line), e = 0.1 (light blue dotted line), e = 0.6 (yellow continuous 
line). For the stellar and planetary physical parameters, and for the inclination of the orbit to the plane of the sky and to the star/planet 
equators we used the standard values of Table [2] We fixed the periastron at uj = 90 deg. 



is due to the overall change of y over one orbital revolution. 
Since the linear speed of the guiding center is va = na, 



P tr = P b + ^ = ^(l + Atr ). 
vq n 



(80) 



with 



Ati 



Aj/_ 
2-7ra' 



(81) 



the approximation of eq. (|79[) and 12 l 



!) 



lae sin /, 



(82) 



it is possible to explicitly compute the effects of several dy- 
namical perturbations on P tI . The averaged time variation 
of it can straightforwardly be computed from eq. |T]) for 
Y — )> y /na, with y conveniently given by eq. (|82[) . A com- 
mon feature of the results that we are going to show below 
is that non-zero long-term variations of P tI occur for all the 
dynamical effects considered; they do not vanish in the limit 
e^0. 



in which Ay is the variation of y per orbit, i.e. integrated 
over one orbital period. By using the first term of order 0(e) 



12 As a consequence, p moves about G in the opposite sense 
with respect to G abou t F on a 2:1 ellipse of period 2n/n 



JMurrav &: Dermotd feoOOl . 
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6.1 The effect of the stellar oblateness 

A straightforward calculation yields the following non-zero 
long-term harmonic effect. 



dP tl \( J D 

dt I 



' R *\ 2 V 
TTHT \-T) X 



64(l-e 2 ) 3 

x [-(2 + e 2 )(-8 + 5e 2 )(l + 3cos2**)~ 



;+ ?e 2 ) sin 2 ** cos 2oj1 



(83) 
Note that eq. (|83[) vanish neither for equatorial (^* = 0) 
nor for polar (<]/* = 90 deg) orbits. 



6.2 The tidal bulges 

It is found that the tidal bulges cause non-vanishing long- 
term seculaa 13 ! variations of P tr . They are 



dP tr I ( tid P) 

dt I 



(—) (-) 



r p s5 3(-128-464e 2 +320e 4 + 195e 6 )fc 2p 
"^ 64(l-e i! )S ' 

(84) 

dPtr I (tid *) \ _ / mp \ , Rt ^5 3(-128-464e 2 +320e 4 + 195<: 6 )k2* 

"dt^l / ~~ V^*/ ^~' 64(l-e 2 ) 6 

(85) 
Note that, also in this case, the ratio of the planet-to-star 
tidal effects is equal just to eq. (|58[) . as in Section ^. 3l for V p 
and Section \5. II for Ai C ci- 



6.3 The effects of general relativity 

The general relativistic gravitoelectric term induces the fol- 
lowing long-term secular variation of P tr . 



dPu 
dt 



(GE) 



He 



(24 + 33e 2 - 7e 4 
4(1 -e 2 ) 2 



(86) 



It is interesting to note tha t eq. (|86[) is quite different from 
eq.(17) bv ljordan fc Bakosl (J2008) yielding an instantaneous 
formula proportional to e. Also the gravitomagnetic field 
causes a net non- vanishing harmonic effect, which is 



dt 



(GM) 



GS* 



i'-i) 



cos *!/* 



(l_ e a)3/a 



(87) 



It vanishes for polar orbits. 

6.4 The third body X 

As expected, in this case the calculations are much more 
cumbersome. To effectively perform them, the further ap- 
proximation 

(1 + ecos/) - « 1 — fcecos/, k = 2, 3, 4 

is adopted. We have 



dPu 
dt 



(X) 



/C x \ 3(l-e 2 ) 2 2J(f x ,n,^,*,) 

^ 16 ' (89) 



where 2J is a complicated function of the angular orbital 
parameters of the planet and of the position of X. As in 

13 Indeed, a and e do not undergo long-term time variations. 



the previous cases, it is not worth explicitly showing it. It 
is found that the long-term harmonic signal of eq. (|89[) does 
not vanish for any particular orbital configurations of both 
the perturbed planet p and the distant perturber X. For 
example, if both lie in the equatorial plane of the hosting 
star, eq. (|89p is not zero. 



6.5 Numerical evaluations 

The standard scenario of Table [2] yields 

= 1 x 1(T 6 , 



dP tI I (GE) 

dt I 



dP tr I ( tid P) 

dt I 



dPtr I OP 

dt I 



dP tT \<.J%) 

dt I 



dPtt I (tid *) 

dt I 



dP tI I (GM) 

dt I 



= 8x10"', 

= 1 x 10~ 7 , 

= 8xl0~ 9 , 

= 4x 10~ 9 , 

= 1 x 10" 10 



(90) 



The effect of a Jupiter-sized third body X at rx = 0.1 — 1 au 
is of the order of 10~ 5 — 10~ 8 . An Earth-type perturber X 
at the same distances would yield variations of Ptr of about 
10 -8 — 10 -11 . Figure \5\ depicts the magnitude of the upper 
bound of < Ptr ) for all the dynamical effects considered, with 
the exception of X, as a function of the semi-major axis for 
different values of the eccentricity. 



6.6 Measurement accuracy 

Concerning the accuracy with whic h time variations of the 
transit period could be detected, iMiralda-Escu dc (2002) 
yields ap ~ lO^ 1 over t — 10 yr. He assumes about 
OTiti = 1000 primary transits and an accuracy in measur- 
ing each primary transit of cri ec i ~ 1 s . A similar evalu- 
ation is given by Ijordan fc Bakosl i|2008l h Note that their 
formula of eq.(17) yields a value for the general relativistic 
gravitoelectric effect which is about 6 orders of magnitude 
sm aller than ours in eq. (|90p ; the pessimistic conclusions 



by Ijordan fc Bakosl 1120081) are actually based just on such 
a sma ll figure. Finally, let us note that iMacieiewski et alj 
(2010) have recently detected transit timing variations in 
WASP-3b (a = 0.03 au, e w 0, P b = 1.8 d) over r = 2 
yr (2009-2010), but they exhibit a harmonic pattern with a 
periodicity of about 124 d and semi-amplitude of 1.4 x 10~ 3 
d, with no discernable secular trends. In fact, such kind 
of cumulatively growing patterns, which may be caused by 
the tidal bulges and the general relativistic gravitoelectric 
field according to eq. (|84|l -eq. ([84} and eq. (|86fl . could not 
have been detect ed since the accuracy wi th which P tr has 
been measured bv lMacieiewski et al.l 1)20101) is of the order of 
2.4Pb = 4 x 10 5 s. Over r = 2 yr this naively translates into a 
potential accuracy in detecting P tr of about 6 x 10 -3 . Look- 
ing at eq. (|89[1 for X, it results that long-period harmonic 
effects are, in principle, possible in view of the interplay 
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Figure 5. Maximum values of the long-term time variations ( Ptr) as a function of a (0.014 au < a < 0.449 au) for different values of 
the eccentricity: e = 0.005 (red dash-dotted line), e = 0.03 (blue dashed line), e = 0.1 (light blue dotted line), e = 0.6 (yellow continuous 
line). For the stellar and planetary physical parameters, and for the inclination of the orbit to the plane of the sky and to the star/planet 
equators we used the standard values of Table [2] We fixed the periastron at uj = 90 deg. 



between the node and the periastron of WASP-3b entering 
2); recall that they slowly change because of the oblateness, 
tidal bulges, general relativity and a X themselves. However, 
the figures for ( Ar) in Section fo.SI valid for plausible values 

of /Cx, do not allow to obtain magnitudes as large aa 14 l 1CP 3 
d. 



7 SUMMARY AND CONCLUSIONS 

We looked at the impact of several classical and general 
relativistic dynamical perturbations on directly observable 



They can be obtained by multiplying { P% r ) by Pj, . 



quantities in transiting exoplanets. We considered the cen- 
trifugal oblateness of both the hosting star and the planet 
due to their rotations, the tidal bulges mutually raised by 
both the star and the planet on each other, and a distant 
third body X as Newtonian effects. Concerning general rel- 
ativity, we took into account both the so called gravitoelec- 
tric, Schwarzschild and the gravitomagnetic, Lense-Thirring 
perturbations. The observables considered are the transit 
duration Atd, the radial velocity V p , the time interval At ec i 
elapsed between primary and secondary eclipses, and the 
time span Ptr between successive primary transits. 

We analytically worked out, in an uniform and straight- 
forward way, the long-term, i.e. averaged over one orbital 
revolution, temporal variations of such quantities caused by 
the aforementioned non-Keplerian features of motion. We 
did not restrict ourselves to circular orbits. 
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For all the dynamical effects considered we released nu- 
merical evaluations of their magnitudes by adopting a typ- 
ical exoplanetary scenario involving a Jupiter-sized planet 
closely (a = 0.04 au) revolving around a main sequence Sun- 
like star along a moderately eccentric orbit (e = 0.07) lying 
in the orbital plane of the parent star (** = 0). The rotation 
of the planet has been assumed synchronized with its orbital 
frequency because of tidal effects, and an edge-on configu- 
ration (i — 90 deg) has been assumed. For each observable 
considered, we gave an order-of-magnitude evaluation of the 
accuracy with which it could be measured over an observa- 
tional time span r = 10 yr. In view of the wide distribution 
of orbital configurations which is expected to characterize 
the numerous planets which will be discovered by the ongo- 
ing space-based Kepler mission, we also performed graphical 
investigations of the dependence of the effects studied on the 
semi-major axes a and the eccentricities e. 

No net changes ( Atd ) in the duration transit occur for 
both the tidal bulges and the Schwarzschild perturbations. 
Instead, the oblateness of the bodies, a distant planet X 
and the Lense-Thirring effect induce long-term, harmonic 
variations of the transit duration which do not vanish even 
for circular orbits. For exactly edge-on orbits (i = 90 deg), 

(Atd) — in all cases. The same occur, independently of 
the inclination i of the orbital plane to the line of sight, for 
equatorial orbits (>]/* = 0), apart from X if it is not coplanar 
with the perturbed body. If both the planet and its perturber 
X lie in the same plane coinciding with the equatorial plane 
of the star, then ( Atd ) = 0. By assuming small departures 
from the equatorial (\&* = 15 deg) and edge-on (i = 87 deg) 
orbital configuration, the magnitude of the upper limits of 
the non-zero oblateness and gravitomagnetic effects is of the 



order of 10 



-10" 



A jovian perturbing body X at 0.1 — 1 



au would cause a rate of ~ 10 — 10 , while for an Earth- 



size X at the same distances we have < At 



io~ 7 -io" 



The expected accuracy in measuring ( Atd ) over r = 10 yr 

is about a k t ^ 10~ 9 - 10~ 8 . 

The radial velocity V p undergoes non-vanishing long- 
term, harmonic variations caused by all the perturbations 
considered only for eccentric orbits; the resulting amplitudes 
are proportional to e, and vanish neither for equatorial nor 
edge-on orbits. The largest effects, of the order of 10~ 7 m 
s~ 2 , are due to the tidal bulge raised on the planet and 
the general relativistic Schwarzschild term. The centrifugal 
oblateness of the planet causes an acceleration of the order 
of 10 -8 m s -2 , while the changes due to the centrifugal 
oblateness and the tidal distortion of the star are of the 
order of 10~ 9 m s~ 2 . The Lense-Thirring effect amounts to 
The perturbation of a distant Jupiter-sized 



body X at 0.1 — 1 au is ( V p 



Kl- 



lO m s , while a 



rocky planet with the mass of the Earth located at the same 



distances would cause ( V, 



io- 



10" 



The 



accuracy with which it is possible to measure secular trends 



in the radial velocity is a^ 



10 



9 m s 2 over r = 10 yr. 



Concerning the time interval At cc i between primary and 
secondary eclipses, it experiences long-term, harmonic vari- 
ations only for eccentric orbits. Its amplitudes are propor- 
tional to e, and are non-vanishing for all the perturbations 
considered. Also in this case, the largest effects, of the order 



of 10 -7 , are due to the tidal bulge raised on the planet and 
the general relativistic gravitoelectric field. The centrifugal 
oblateness of the planet causes an effect of the order of 
10~ 8 , while the variations induced by the centrifugal oblate- 
ness and the tidal distortion of the star are of the order of 
10~ 9 . The effect due to the general relativistic gravitomag- 
netic field is as large as about 10 -11 . The perturbation of a 
jovian-type body X at 0.1 — 1 au is /At ec i) ~ 10~ 5 — 10~ 8 , 
while a terrestrial planet at the same distances would in- 



duce (At 



10" 



io- 



The expected accuracy in 



measuring ( At ec i ) over r = 10 yr is about a^ t ~ 10 

All the dynamical effects considered causes long-term 
variations of the transit period P tr . They do not vanish 
for circular orbits. While the signals due to the stellar and 
planetary oblateness, the gravitomagnetic field and a third 
body X are harmonic, those induced by the tidal bulges 
and the gravitoelectric field are secular rates. The largest 
effect, of the order of 10 -6 , is the trend due to general rel- 
ativity (Schwarzschild). The rate due to the planetary tidal 
bulge is 8 x 10~ 7 , while that due to the star's tidal distor- 
tion is 2 orders of magnitude smaller. The amplitudes of 
the signals caused by the planetary and the stellar oblate- 
ness are 1 x 10~ 7 and 8 x 10~ 9 , respectively. The gravit- 
omagnetic effect has an amplitude of the order of 10~ . 
The action of a distant Jupiter-like planet X at 0.1 — 1 au 
P tl > « 10~ 5 — 10 -8 , while an Earth-type planet at the 



same distances would cause < Ptr 
be possible to measure long-term changes in P tI with an ac- 
curacy of about 10~ 13 over r = 10 yr. 

Our results may be useful in focussing the research on 
some specific exoplanets or classes of exoplanets, particu- 
larly tailored to measure a given dynamical effect of in- 
terest. Moreover, they should allow to better evaluate the 
systematic uncertainty in the determination of a given stel- 
lar/planetary physical parameter induced by the other per- 
turbations which are to be considered as sources of bias. 
Moreover, our approach can, in principle, be extended also 
to other scenarios and dynamical features of motion involv- 
ing, e.g., natural and artificial bodies in our solar system in 
order to suitably design new tests of general relativity and 
modified models of gravity. The latter issues, along with 
a more detailed quantitative analysis of different, specific 
exoplanetary scenarios, may be the subject of further re- 
searches. 
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